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ABSTRACT 

The development of probabilistic structural analysis methods is a 
major part of the SSME Structural Durability Program and consists of 
three program elements: (1) composite load spectra, (2) probabilistic 

finite element structural analysis, and (3) probabilistic structural 
analysis applications. Recent progress includes: (1) the effects of 

the uncertainties of several factors on the HPFP blade temperature 
pressure and torque, (2) the evaluation of the cumulative distirbution 
function of structural response variables based on assumed uncertain- 
ties on primitive structural variables, and (3) evaluation of the 
failure probability. collectively, the results obtained demonstrate 
that the structural durability of critical SSME components can be 

probabilistically evaluated. 


INTRODUCTION 


It is becoming increasingly evident that deterministic structural 
analysis methods will not be sufficient to properly design critica 
structural components for upgraded Space shuttle Main Engines ( 
Structural components in the SSME are subjected to a variety of Com- 
dex severe cyclic and transient loading conditions including high 
temperatures and high temperature gradients. Most of these are quan- 
tifiable only as best engineering estimates. These complex loading 
conditions subject the material to coupled nonlinear behavior which 
depends on stress, temperature, and time. Coupl ed nonl inear mate ri 
behavior is nonuniform, is very difficult to determine experimentally, 
and perhaps impossible to describe deterministically. In addition, 
test data on critical SSME structural components are relatively small* 
Fabrication tolerances on these components, which in essence are small 
thickness variations, can have significant effects on the c °™P° n ^ 
structural response. Fabrication tolerances by their very nature are 
statistical. Furthermore, the attachment of the components to t e 
structural system generally differs by some indeterminant degree from 
that which was assumed for designing the component. In summary, a 
four fundamental aspects of: (1) loading conditions, (2) nvaterial 

behavior, (3) geometric configuration, and (4) supports - on which 
structural analyses are based, are of a statistical nature. One direct 
way to formally account for all these statistical aspects is to develop 
probabilistic structural analysis methods where all participating 
variables are described by appropriate probabilistic functions. 

NASA Lewis Research center is currently development probabilistic 
structural analysis methods for select SSME structural components under 
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the SSME Structural Durability Program. Briefly, the deterministic, 
three-dimensional, inelastic analysis methodology developed under the 
Hot Section Technology (HOST and B&T Base Programs) is being augmented 
to accommodate the complex probabilistic loading spectra, the thermo- 
viscoplastic material behavior, and the material degradation associated 
with the environment of space propulsion system structural components 
representative of the SSMF, such as turbine blades, transfer duct, and 
liquid-oxygen posto, Fig. 1. 


The development of probabilistic structural analysis methodology 
consists of the following program elements: (1) composite load 

spectra, (2) probabilistic structural analysis methods, and (3) proba- 
bilistic structural analysis application (ref. 1). The program main 
goal is to develop the methodology to address the problem depicted 
schematically in figure 2. Past progress of tasks in specific elements 
of the program are described in papers presented in conferences 
(ref. 2, 3, and 4). 


Pecent activities focused on extending the methodology to include 
the combined uncertainties in several factors on the structural 
response. An executive summary of this progress is shown in figure 
3. The objective of the present paper is to briefly describe progress 
in three program elements: composite load spectra, probabilistic 

finite element structural analysis, and strength degradation. Progress 
is described in terms of fundamental concepts, computer codes, and 
re pr esen ta t ive results. 


COMPOSITE LOAD SPECTPA 

The fundamental assumption for developing composite load spectra 
is that each individual load condition is the probabilistic time 
synthesis of four primitive parts: (1) steady state, (2) periodic, 

(3) random, and (4) spike. Each of these parts, except random, is 
described by a deterministic portion and a probabilistic perturbation 
about this deterministic portion as depicted schematically in figure 
4. One justification for describing each loading condition in terms 
of primitive variables is that: experts over the years have developed 

good judgments of the ranges of the perturbations about nominal 
(deterministic) conditions. The objective of the Composite Load 
Spectra program is to formalize the fundamental assumption in a 
computer code using: (1) available data from various rocket engines, 

(2) probability theory, and (3) a dedicated expert system. 

A schematic diagram of the composite Load Spectra (CLS) Computer 
code is shown in figure 5. Representative results obtained for the 
perturbations of different engine factors on the high pressure turbo- 
pump blade are shown in figures 6 and 7. Figure 6 is the nominal 
temperatures, while figure 7 indicates the temperature changes due to 
hot gas seal geometry and respective perturbations indicated in the 
figure caption. For example, the greater temperature change due to Gas 
Seal Geometry (fig. 7) of 0.06 is 53.3 °F. The combined contributions 
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of this and other factors (not shown here) is 87.6 °F which is in addi- 
tion to the greatest nominal temperature of 1860 °F. Although at 

first glance a change of 87.6 °F may seem insignificant, this is not 

the case because at these high temperature small temperature changes 
have dramatic effects on the material structural durability and atten- 
dant cooling requirements. 

Another representative example is shown in figure 8, where the 
effect of comparable factors on the torque of the High-Pressure- 
Ox i dizer -Turbopump ( HPOTP ) are plotted as bounds versus time. These 
bounds are substantial at some times and relatively close at others. 
Similar plots can be obtained for pressures or any other loading con 
dition. The current CLS code permits the simultaneous perturbation of 

47 different factors for each different load condition. 


PFOBABIL I STIC FINITE ELEMENT STFUCTUFAL ANALYSIS 


The fundamental assumption for developing probabilistic finite 
element structural analysis (PFESAN) is that "the uncertainties in 
each primitive structural variable can be described by assumed corres- 
ponding probabilistic distributions." Primitive structural variables 
are those which are used to describe a structure such as: (1 

ness, (2) strength, (3) thickness and tolerance, (4) spatial location, 
(5) attachment, (6) various nonlinear dependencies (temperature, 
stress, time, etc.). A schematic of the probabilistic distributions 
for some primitive variables is shown in figure 9. subsequently, the 
uncertainties in the load conditions (described by the composite load 
spectra) and the uncertainties in the primitive structural variables 
are computationally simulated by performing multiple finite element 
structural analysis to determine the probabilistic structural response 
of a specified SSME structural component. The structural response is 
generally described in terms of displacement, frequencies, buckling 
loads, and structural fracture toughness. The integration is ill us 
trated schematically in figure 10. 

It is instructive to compare component development by the tradi- 
tional engineering approach and component evaluation using PFESAN. 
The par all ism is summarized in Table l. The former approach relies on 
physical experimental and requires that the material, fabrication pro 
cess, and test methods are already available. The latter approach 1 
entirely computational and requires the integration of . availa le 
structural analysis methods with available probability theory. The 
former approach has the advantage of demonstrating a specific tec 
nology while the latter has the advantage of assessing undeveloped but 
with high payoff potential candidate technologies. In addition 
former approach requires a large number of experiments while the 
latter can be verified with strategically selected few. 


PFESAN has been formalized 
NESSUS (Nonlinear Evaluation of 
NESSUS is driven by a dedicated 


and integrated into a computer code 
Stochastic Structures Under Stress), 
expert system. A schematic diagram of 
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NESS US is shown in figure 11. The user interacts with NESSUS through 
a dedicated expert system schematically shown in figure 12. Represen- 
tative results obtained using NESSUS are shown in figure 13. The dis- 
tributions assumed for the primitive variables listed in in the table. 
Both the individual and the combined effects of the primitive variables 
on the combined stress (von Mises) are shown in the figure in terms of 
cumulative distribution functions (CDF). The information generated 
during the PSESAN can be used to establish bounds on the CDF. A sample 
result is shown in figure 14 for one blade location. 

The curves in figure 14 can be used in a number of ways. Two of 
them are: (1) all the blades tested in the assumed conditions will 

have a mean combined stress between 57.9 and 62.1 ksi 90 percent of 
the time; (2)the mean combined stress in all the blades tested (under 
the assumed conditions) will range from 42 ksi to 83 ksi. This indi- 
cates that a wide scatter in the mean combined stress is probable. 

Two implications follows; (1) assessing the durability life of the 

blades using only material uncer ta in t i es will not be sufficient, and 

(2) obtaining wide scatter in measured stress/strain magnitudes does 

not indicate test procedure difficulties. 

It is noted that all the NESSUS results presented herein were 
obtained using 50 simulations for each case studied. These simulations 
are relatively small compared to direct Monte Carlo simulation which 
will normally require 1000 and greater simulations. The reduced but 

with comparable accuracy simulations is a NESSUS feature which uses 
the fast probability integration method to select subsequent simula- 
tions in a self-adaptive manner. 

PROBABILISTIC ANALYSIS FOR STRENGTH DEGRADATION 

The fundamental assumption for developing probabilistic analysis 
methods for strength degradation is that "the uncertainties in primi- 
tive variables of the strength degradation model can be described by 
assumed probabilistic distributions." Two different models were 
selected to demonstrate the concept. The models express the number of 
mechanical load cycles to failure. One of the models is based on 
linear elastic fracture mechanics and the other on a strength degrada- 
tion model recently studies at Lewis. The models with their respective 
primitive variables are summarized in table 2. 

Both of these models were used to predict the number of cycles to 
failure in a material used in SSME components. The input for the 
fracture mechanics model is summarized in table 3. The CDF obtained 
from this input is shown in figure 15. The input for the strength 
degradation model is summarized in Table 4 and the corresponding CD F 
is shown in figure 16. Both CDF exhibit wide ranges for the probable 
number of cycles to failure. The linear fracture mechanics model shows 
a mean of 10 ,000 cycles while the strength degradation model shows a 
mean of 10 million cycles. Based on this comparison the linear frac- 
ture mechanics model penalizes the material by three decades. It is 
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important to note the differences between the two models: (l) he 

linear fracture mechanics model assumes the existence of ^ crack-1 ike 

and then evaluates the number of cycles required to grow this 

defect to a critical size for imminent rapid propagation to fracture. 

f 2 The strength de gC adation model does not presuppose the existence 

of defect^ an’d, therefore 

gallon resulting in greater "“^ t ° v VC c Ta b les while the strength 
mechanics model has > P the greater the number of primitive 

d egra b tr;n ha t S he 3 ;odeTrmo?e inclusiv e the representation of the 
variables in the m strength degradation model will be more 

PhY uJate in (4) ^eYinear fracture mechanics model requires determina- 
accurate. (4) 1 specialty and often complex test methods 

whil e° the" 'strength degradation model uses available room temperature 

material properties. 


irrespective of the model used, the important 
the uncertainties in fatigue cycles to failure 
probabil istically . 


conclusion is that 
can be evaluated 


CONCLUSIONS 


The development or probaMl 1st Si 

Ifllfts Tt^rTSinSS*™ “.««1 incurs on blade temperatures. 

and hnraue (2) the evaluation of the cumulative distribu 
5S“^i'tlS .f structural rescue variables based on assumed I uncer- 
tainties in the structural . ^ ! J s a “„ t “ ' teri o£ cumulative 

™e mechanics nn^n -ength 

degradation . model. *^ *"““^^*(0,8) for the probabilistic 

an.lv.ts of select SEME structural co»pc >"«£, and (3. . Uf - ^ 
ity code for the assessment o ^‘“^^ents. collectively, the 

^ItsTbwTed to date demonstrate that the structural durability of 
isME criutal components can be evaluated using the methodology devel- 
oped under the SSME Structural Durability Program. 
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TABLE 1. 


SOME HELPFUL PARALLELS 


CaiPPHENT DEYELPmENT 

PARAMETER IDENTIFICATION 
CHARACTER CAT ION - DATABASE 
COMPONENT FABRICATION 
COMPONENT TESTING 
COMPONENT TEST RESULTS DATABASE 
STATISTICAL INFERENCES 
RELIABE nr /CONFIDENCE LEVEL 


P5AM COMPONENT EVALUATION 

PARAMETER IDENTIFICATION 
PROBABILISTIC DISTRIBUTIONS 
COMPONENT MECHANISTIC MODEL 
COMPONENT ANALYSIS 
COMPONENT ANALYSIS RESULTS OAT/CASE 
STATISTICAL INFERENCE 
RELIABILnY/CONFIDENCE LEVEL 


BKSIttL 


BOTTOM LINE 

COMPUTATIO NAL 


TABLE 2. 

PROBABILISTIC ANALYSIS FOR STRENGTH DEGRADATION MOOELS 

MOOELS SELECTED FOR STUDYi 

1. FATIGUE CRACK GROWTH MODEL (PARIS EQUATION) 

r -•*>/•» 


NW 


a, 


cy w n m/i act ^ 

WHERE C. M. ATT AM) ARE RANOOM VARI/WLES 


2. STRENGTH REDUCTION MODEL (CHAMIS EQUATION) 


' *>/ 2 * I 


- ♦ I 





WHERE S q , T r . T, T Q , $ r . CT . C£ . 

AN) WE RANOOM VWIWLES 


I 

j 
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TABLE 5. RANDOM? INPUT (FATIGUE CRACK GROWTH MODCD 


VARIABLE 

VARIABLE 

TYPE 

DISTRIBUTION 

TYPE 

MEAN 

STANDARD 

DEVIATION 


C (MATERIAL 
PROPERTY) 

RANDOM 

LOGNORMAL 

8.66 x lO'IO IN/CYCLE 

0.866 x lO'IO 

(101) 

N (MATERIAL 
PROPERTY) 

RANDOM 

NORMAL 

6.M8 x 10'? IN/CYCLE/, 
KPSI INW? 

3.?HxlO-J 

(St) 

A^( ALTERNATING 
STRESS) 

RANDOM 

LOGNORMAL 

90 KPSI 

9 

not) 

A, (INITIAL 
CRACK SIZE) 

RANDOM 

LOGNORMAL 

118 x lO'R IN 

17.7 x lO'R 

MSI) 

• r (FINAL _ 
(CRACK SIZE) 

DETERMINISTIC 

N/A 

7.87 X 10'? IN 

N/A 


T (COMPONENT/ 
CRACK SHAPE PAR.] 

DETERMINISTIC 

1 

N/A 

1.0 

N/A 



TABLE, A. RANOCHJ AND RANDOMM INPUT (STRENGTH REDUCTION MODEL) 


VARIABLE 

VARIABLE 

TYPE 

DISTRIBUTION 

TYPE 

MEAN 

STAND ARO 
0EV1AI10N 

T r (MELTING TEMPERATURE) 

RANDOM 

NORMAL 

2732*F 

82. 

(31) 

Sr (ULTIMATE TENSILE STRENGTH) 

RANDOM 

LOGNORMAL 

130 KPSI 

6.S 

(SI) 

Nnr (LOG OF FINAL CYCLE) 

RANDOM 

LOGNORMAL 

6 

0.8 

(101) 

T 0 (REFERENCE TEMPERATURE) 

RANDOM 

NORMAL 

68 *F 

2.0 

(31) 

</£ (RESIDUAL COMPRESSIVE STRESS) 

RANDOM 

lognormal 

-?.9 KPSI 

-0. IMS 

(SI) 

Ifao (LOG OF REFERENCE CYCLE) 

RANDOM 

LOGNORMAL 

7 

0.7 

(101) 

S* (REFERENCE FATIGUE STRENGTH) 

RANDOM 

LOGNORMAL 

72.6 KPSI 

3.6 

(SI) 

T (CURRENT TEMPERATURE) 

RANDOM 

NORMAL 

1562*F 

M6.7 

(31) 

O' (CURRENT MEAN STRESS) 

RANDOM 

lognormal 

21.8 KPSI 

1.1 

(SI) 

S (CURRENT FATIGUE STRENGTH) 

RANDOM 

LOGNORMAL 

36.3 KPSI 

1.8 

(SI) 

N (TEMPERATURE EXPONENT) 

RANDOM 

NORMAL 

O.S 

0.15 

(31) 

n (STRESS EXPONENT) 

RANDOM 

NORMAL 

o.s 

0.15 

(31) 

0 (CYCLE EXPONENT) 

RANOOM 

NORMAL 

O.S 

0.1S 

(31) 
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FIGURE ! . SELECT SSME COMPONENTS 
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FIGURE 2. PROBABILISTIC STRUCTURAL ANALYSIS DEFINITION 
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FIGURE 3. PROBABILISTIC STRUCTURAL ANALYSIS METHODS (PSAM) 
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FIGURE <t. 

EACH LOAD CONDITION IS THE PROBABILISTIC TIME SYNTHESIS OF FOUR 
PRIMITIVE PARTS: COMPOSITE LOAD SPECTRA 
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FIGURE 5. 

COMPOSITE LOAD SPECTRA SIMULATION USING EXPERT SYSTEMS 

10EXPT DESIGN 

COMPOSITE LOAD SPECTRA EXPERT S VS TEN LOAD EXPERT SYSTEM DESIGN PHILOSOPHY 
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FIGURE 6. COMPOSITE LOAD SPECTRA PREDICTED STEADY STATE TEMPERATURES 
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FIGURE 7 TEff>ERATURE PERTURBATION DATA DUE TO 
r '• 0ELTP temperatures 

ONE STANDARD OEUIATION OF .060 
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FIGURE 8. TRANSIENT MODEL HPTOTP TORQUE 




FIGURE 9. 


PROBABILISTIC STRUCTURAL ANALYSIS: THE UNCERTAINTIES IN EACH 
STRUCTURAL PRIMITIVE VARIABLE CAN BE DESCRIBED BY ASSUMED 
CORRESPONDING PROBABILISTIC DISTRIBUTIONS 
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FIGURE 10. 

THE UNCERTAINTIES IN THE LOAD CONDITIONS AND STRUCTURAL 
PRIMITIVE VARIABLES ARE INTEGRATED BY PERFORMING MULTIPLE 
FINITE ELEMENT ANALYSES 
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FIGURE 11. 

PROBABILISTIC FINITE ELEMENT COMPUTER CODE NESSUS 
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FIGURE 13. PROBABILISTIC STRUCTURAL ANALYSIS USING NESSUS 
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FIGURE 16. STRENGTH REDUCTION MODEL 
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